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ABSTRACT 


The modulated beam of the free electron laser when 
passed through a thin target should produce strong 
transition radiation. The measurement of the transition 
radiation is a direct measurement of the electron beam 
modulation in the free electron laser. A transition 
radiation experiment using the Stanford MKIII Infrared Free 
Electron Laser (IRFEL) has been proposed. The analysis has 
centered on TRANSPORT, a computer program used for designing 
charged particle beam systems. The MKIII IRFEL wiggler exit 
bending magnet system was modeled using TRANSPORT. Analysis 
reveals that the transverse emittance and momentum spread 
characteristics will cause the modulated beam to demodulate 
along the path of the central trajectory. A detector 
location 10 millimeters downstream of the first bending 
magnet is found suitable for the measurement of the electron 
beam modulation. For this case the thin foil must be 
rotated approximately 11.25 degrees about the y-axis to 
yield a minimum effective picobunch extent as seen by the 
Eady mero: Alternate FEL systems are suggested for 


investigation. 
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I. INTRODUCTION 


A. OVERVIEW 


The electron beam from a free electron laser (FEL) 
offers the possibility of observing and using transition 


radiation (TR) at optical frequencies. After interacting 


with the wiggler and radiation fields, the electron beam 


becomes modulated in the axial direction at the lasing 


frequency [{Ref. 1]. The axially modulated bunches 


(picobunches) could be passed through a thin foil producing 


the surface effect TR as illustrated in Figure 1-1. 


ie viewport 
backward/ 


specular TR 
— 


w forward TR 


picobunches 





Figure 1-1 Transition Radiation (TR) 


The axially modulated beam picobunches impact a thin foil of 
aluminum or mylar yielding forward and backward/specular TR 
[Ref. 2]. Both forms of radiation should be enhanced by 
coherent emission from all charges in the picobunch. The 
backward/specular TR is redirected to a viewing port/optical 
detector by a mirror. The resulting radiation could be used 
to diagnose the electron beam modulation [Ref. 1]. 

The TR power radiated per picobunch fie found to be 
proportional to the charge of the picobunch squared. Fora 
picobunch of approximately 10° electrons, the radiation 
intensity can easily be detected and measured with an 
Optical “detececor: The axially modulated electron beam of 


the IRFEL, upon impact with the thin foil or multifo1l) yee 


generate coherent optical transition radiation. [Ref. 2] 


B. MODULATION OF THE ELECTRON BEAM 

The longitudinal bunching of the ,eteerron beam by the 
undulator of the free electron laser (FEL) is a result of 
the force acting on the individual particles by the electric 
field of the superimposed laser radiation. This effect ona 
Single electron is examined in Figure 1-2 (Ref. 3]. This 
figure illustrates the fields and trajectories within the 
undulator of a FEL: an alternating magnetic field, a co- 
propagating laser electric field and an electron along a 
trajectory. The electric field of the laser is zero at 
position A. The electron experiences no force other than 


that of the undulator's magnetic field, which will bend it 








electric field of 
laser radiation 


undulator — Ve Crime 
magnetic field (may ecrony 





Figure 1-2 Longitudinal Bunching of the Electron Beam 


in the horizontal plane. Now consider the electron a time 
later at position B. The electron has traveled one-quarter 
of an undulator period and now experiences a maximum in the 
laser electric field. The force on the electron at point B 
is opposite the direction of the laser electric field at 
that point and is labeled F,. The electron will therefore 
experience a decelerating force. Whether an electron is 
accelerated or decelerated by the laser electric field 
depends on the phase of the laser electric field with 
respect to the electron oscillation. Some electrons are 
accelerated, while those half a laser wavelength ahead or 
behind are decelerated. For example, an electron at C would 


be accelerated. As the accelerated electrons catch up to 


those which have been decelerated, they form into bunches at 
the laser wavelength. In other words, the electron beam 
becomes axially modulated at the laser wavelength of the FEL 
radlatlonseeirRer.- 3a 

Figure 1-3 [Ref. 4] depicts approximately 10° electrons 
being bunched within the optical wavelength of the FEL 


radiation while within the undulator. 


Le 
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Figure 1-3 Bunching of Electrons within the Wavelength 
of Laser Radiation 


As the electrons oscillate they spontaneously emit "magnetic 
bv ene eeeanmunse peaked in the forward direction thus 
generating the laser electric field. When the laser 
frequency 1s nearly resonant with the electron oscillations, 
Significant amounts of energy can be exchanged between the 
electrons and the laser field causing the bunching effect of 
the electron bean. The condition for this resonance is 


given by Equation 1-1 (MKS) [Ref. 5]: 


where: 


y = energy of resonant electron divided by rest 
mass energy, 


My = wiggler (undulator) magnetic period, 
X = laser wavelength (operating wavelength), 
e = electron charge, 
Bems = rms magnetic induction of the wiggler, 
ee — velocity of light in vacuo. 


The operating (lasing) frequency of a FEL is determined 


By Equation 1-2 (MKS), 
d eee ama KO), ae” 


where: 


K = the dimensionless measure of the wiggler 
strength, 


The operating parameters of the Stanford MKIII IRFEL are 
shown in Table 1-1 [Ref. 6]. 

The wiggler period of the IRFEL is of the "hybrid" type 
employing both SmCo permanent magnets and vanadium permendur 


steel poles. The magnets are oriented to provide a nearly 


TABLE Fi 


MKIII IRFEL PARAMETERS 


WIGGLER 

magnet period Aw Qed cm 
number of periods N 47.5 

minimum gap Smin 655 mm 
Wiggler strength K 14 

peak magnetic field Bo 0.702 7 
RESONATOR 

length L l7 Soe em 
Rayleigh range ZR 73 cm 
E-BEAM 

energy Y 85 

energy E 43 MeV 
energy spread Y/Y O05 % 
peak current I 20 Amps 
emittance En 7 T-mm mrad 
macropulse length Uv 3 usec 


Sinusoidal vertical field By, using a Halbach undulator 
design. (Ref. 6] Tne lasing frequency can be modified by 
slightly changing the energy of the beam or by changing the 
K value of the wiggler. In practice the K parameter is 
modified by adjusting the gap between the permanent magnets 
in the wiggler cavity. The resonator length must then be 


modified to account for the new wavelength. Assuming the 


lasing wavelength of 3.0 microns, the IRFEL will produce a 


modulated beam as shown in Figure 1-4. 





Figure 1-4 Disks of Charge (picobunches) 


The rectangles indicate disks of charge of thickness 1 
micron separated by 3 microns. This periodic bunched pulse 
train will generate transition radiation of wavelength 3 
microns and also odd harmonics when passing through a thin 


foil such as mylar or aluminum. The current vs time profile 


is shown in Figure 1-5. 





Figure 1-5 Current pulses 


C. WHAT IS TRANSITION RADIATION 
Transition radiation (TR) is the electromagnetic radiation 
which is emitted when a uniformly moving charged particle 
traverses from one medium into another with a different 
dielectric constant. More generally the effect will occur 
in the presence of a heterogeneity in a medium. [Ref.7] 

TR may be considered to originate from the change in the 
fields as an electron crosses a boundary. The magnetic 
field remains essentially constant but the electric field 
changes because of the change in the dielectric constant. A 
third field, radiation, is required to meet the boundary 
conditions. 

TR can also be explained in terms of a charge and its 
image. While an electron moves uniformly in vacuo it does 
not radiate; the field in vacuo is equal to the field of the 
electron and its image moving toward it. When the electron 
crosses the boundary surface both the electron and its image 
cease to exist from the point of view of the field. The 
crossing of the boundary must give rise to exactly the same 
radiation as that which is due to a sudden stop at the same 
point of the electron and its image. [Ref. 7] 

The free electron laser will longitudinally modulate the 
electron beam in the undulator/wiggler yielding disks of 
charge of small thickness (lu). Due to the low efficiency 
of OTR, 10°© electrons per picobunch will yield approximately 
104 photons per picobunch upon foil impact. 


For a relativistic beam of electrons (8g- 1), and a TR 


foil inclined at 45 degrees with respect to the beam axis 


(Figure 1-1), the backward/specular reflection intensity for 


a vacuo/dielectric interface is given by the relation [(Ref. 


Zs 
d2t/dwaQ = F(¥,w) (e4/ Tc) [62/(y7* + 94), (1-4) 
where: 
2 
feaeane enn)” - € cos (Wt6)) 
F(y,w) = I a 2 
[e-sin (+6) ] + € cos (Wté) 
F = Fresnel coefficient for reflection of light 
with the E-vector in the plane of incidence, 
W = angle of specular reflection (45 degrees), 


dielectric constant of foil material. 


m 
it 


The intensity of this radiation should be sufficient for 
the detection and observance of TR from a single foil/multi- 
electron interface. Measurement of the radiation will be 
accomplished via a detector of the phototube, TV camera, or 
reticon type. 

The characteristics of TR that will apply to the 
proposed experiment are: 

(1) The angular distribution of TR is such that for 
relativistic particles the emission angle of the 
photons with respect to the electron beam is 


appproximately equal to l1/y. 


(2) The intensity of TR is a strong function of the 
energy of the particle producing it. 


(3) 


(4) 


(5) 


(6) 


(7) 


The radiation is polarized with the electric field in 
the plane containing the particle velocity vector and 
the wave propagation vector. 


In the case of a periodic medium (i.e., multifoil) an 
rv phase difference is the condition for constructive 
interference or resonance between radiation produced 
at successive interfaces; r 1S a positive integer. 
For m boundaries the intensity varies as m°@. 


If the medium is a transparent dielectric then 
e'=n“(w), where n is the refractive index and e! is 
the dielectric constant. 


For a relativistic electron, TR will arise at the 
surface of the interface. [Ref. 7] 


The power radiated from a transition foil is 
proportional to the charge squared [{Ref. 2]. 
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II. PROPOSED EXPERIMENT USING MKIIIT IRFEL 


A. EXPERIMENT PROPOSAL 

The purpose of this experiment will be to measure the 
electron beam modulation of the FEL using a TR foil. 

The MKIII IRFEL, located at the High Energy Physics Lab 
at Stanford University may lend itself to an experiment to 
measure the axial modulation of the electron beam of the 
FEL. The "bunched" beam downstream of the wiggler magnet is 
normally bent away from the downstream cavity mirror and 
directed to the beam dump. If the beam were to retain its 
axial modulation, upon the passage through a thin foil it 
would generate coherent transition radiation. The thin foil 
will hereafter be referred to as the detector. The detector 
must be conveniently located downstream of the wiggler. 

1. System Description 

The MKIII IRFEL is shown in Figures 2-1, 2-2, 2-3, 
and 2-4. Figure 2-1 is a photo of the accelerator section 
of the IRFEL. Figure 2-2 is a representation of the IRFEL 
downstream of the accelerator. The electron beam produced 
by the accelerator is directed around the upstream optical 
cavity mirror and through the wiggler magnet. The beam is 
axially modulated within the wiggler and is then directed 


around the downstream movable mirror to the beam dump. An 


al 





Figure 2-1 MKIII IRFEL Accelerator 
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external photo of the MKIII IRFEL wiggler magnet is shown in 


Figure 2-3. 





Figure 2-3 MKIII IRFEL Wiggler 


Figure 2-4 is a photo of the experimental beamline apparatus 
downstream of the wiggler cavity. The electron beam will 
proceed from left to right. The figure shows the waveguide 
type piping, flanges, bending magnet #2, bending magnet #3, 
and a four way cross. Bending magnet #1 is upstream to the 
left of the flange and bellows (not shown in the photo). 


The detector could easily be placed on the waveguide piping 


or at the four way cross location. 
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Figure 2-4 MKIII IRFEL Wiggler Exit 


2. Possible Detector Locations 
Figure 2-5 is a top view representation depicting 
possible locations for the detector. The electron bunches 
exiting from the wiggler are deflected by the first bending 
magnet in order to prevent them from striking the resonator 
mirror. 
3. Hardware Design at Detector Locations 
The hardware design needed for a detector between 
locations B and C is shown in Figure 2-6. Placement of the 
detector between points B and C would require modifying the 
waveguide piping, bellows, and flange. Location E is the 


alternative location of the detector and would require the 
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Figure 2-5 Detector Locations 
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Figure 2-6 Detector Placement on Waveguide 
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modification of Figure 2-7A. The linear motion feedthrough 
needs to be connected to a sliding sleeve device in order to 
position the TR foil and mirror. A suitable sliding sleeve 
device is shown in Figure 2-7B. In general, it would be 
easier to modify the IRFEL by placing the detector at 


location E. 
i. ANALYSIS OF ELECTRON BEAM TRANSPORT EQUATIONS 


Note: Readers not familiar with the properties of beams and 
the TRANSPORT development of system equations should read 
Appendices A and B before proceeding with this section. 

The development of the electron beam transport equations 
for locations A through E along the central trajectory is 
found in Appendix B. It is desired to find the best 
location in the beamline to place the detector in order to 
generate coherent transition radiation. Of Spareileular 
concern is the longitudinal extent of the picobunches 
(Figure 1-3) as they traverse the beamline. In the ideal 
case the picobunches would retain their axial modulation 
traversing the beamline. This is the case of a non- 
divergent, monoenergetic beam. In the real world, electron 
beams are non-monoenergetic and are divergent. Beams of 
this nature are specified in terms of their phase space. In 
order to observe coherent transition radiation, the 
picobunches must not "debunch" significantly upon traversing 
the beamline. "Significant" debunching for this experiment 


is on the order of a few (>3) microns. If debunching 


ay 
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Figure 2-7A Detector Placement at Location E 
(All Dimensions in Inches) 
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Figure 2-7B Sliding Sleeve Device 
(All Dimensions in Inches) 
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exceeds 3 microns the coherent TR effect will not be 
observed. A simple representation of the beam traversing 
the system is desired. Figure 2-8 is a representation of 
the paths a single electron located at "a" may take in 
traversing the system. System parameters (obtained from 
blueprint) and definitions applicable to the representation 


are shown below the figure. 


Lo 


central 
trajectory 





Figure 2-8 Trajectory of a Single Electron 
Xo = the initial horizontal displacement of the 
electron with respect to the central trajectory. 


89 = the divergence of the electron with respect to 
the assumed central trajectory. 


Sg = the fractional momentum deviation of the 
electron from the assumed central trajectory. 
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The central trajectory parameters are: 


Ly = drift length 1 S27 os nih 
Epo sect itbetengeh 2 2550 7mm 
L = effective pathlength through Beao Sem 
bending magnet 
a = bending angle 11.25 
P09 = bending radius of magnet 195.30 mm 
Pee nical Cdmiehong th 0.0 mm 
E = energy of electron 40 Mev 
fierce bh. = 0; the electron starts at the magnet pole face) 


1. Buskirk System Equations 


The system equations for this model were derived 
by Dr Fred Buskirk. These equations assume the 
Mieetal drift length,. Lo. ands the effective 
pathlength through the bending magnet, L, equal 


zero. The derived equations are: 


ax = (L1/cos#o. fa) 9 Guat (Lysina/cos4@a) 6, ee ay (2-1) 
ds = -(L,sina/cos2a)@, - (Lytan20)6,, (2-2) 
where: 
ax = the change in the horizontal distance of the 
electron with respect to the central trajectory, 
and 


ds = the difference in the pathlength between the 
arbitrary ray and the central trajectory. 


2 


The electron is to traverse the system from point "a" to the 
xy plane of location E. Assuming the fringing fields are 
negligable the electron will bend in the xz plane only. 


Trajectory 1 £(T1) is for the case of the 


monoenergetic (6, = 0) non-divergent (9, = 0) electron 
initially located at x,. Its trajectory coincides with the 
central trajectory as expected. Trajectory 2 (T2) is te 


case of a non-monoenergetic non-divergent electron initially 
located at xX_,.- Since this electron possesses greater 
momentum than the central momentum it will not be bent as 
much by the bending magnets. The pathlength between the 
bending magnets will be shorter thanw@iy If two electrons 
having equal velocities were located at "a" start downstream 
on trajectories Tl and T2, the electron on T2 would have a 
shorter pathlength through the system. At the time when the 
electron on Tl reached the xz plane of location E, the 
electron on T2 would have already passed throught the system 
by a distance equal to the pathlength difference. T2 ie 
also slightly greater than T1 in the x direction. 
Trajectory 3 (T3) is the case of a non-monoenergetic, 
divergent electron initially located at x,. As this 
electron traverses the system on T3, the pathlength between 
the bending magnets will be different from the central 
trajectory and the electron would end up at a different 
location downstream. The three electron trajectories 


illustrate the change in pathlength due to the effects of 
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divergence and percent momentum change. This is analogous 
to the change in the longitudinal extent of a picobunch as 
it traverses the system. 


2. TRANSPORT Equations 
The TRANSPORT equations for location E are derived 


in Appendix B. These equations are: 


Xp(t) = (2p,_5tan a+ L,/cos? a + Lo) 8% 

+ (-2P,4 + 2Q/cos& + LyjSin 0/cos“a)6,, Cet 23) 
Ip(t) = (-2p59 + 2Q/cos4 +L, sin 0/cos?0)6, 

+ (-2p5,a+ 2Potana + Lytan*a)s4 + lo. (M-15) 


It can be shown that the TRANSPORT equations reduce to the 
Buskirk equations letting the effective pathlength through 
the bending magnet, L, and the initial one-half longitudinal 
extent of the picobunch, 1,, equal zero. The variables x 
and 1 correspond to dx and ds respectively in Equations 2-1 
mee 2-2. TRANSPORT analysis requires the phase space 
representation of the input bean. For a normalized 
horizontal beam emittance of 8 mm-mrad (revised MKIII IRFEL 
parameter) with a beam energy of 40 MeV, assume the 2-D 
phase space of Figure 2-9. The points a-d corresponding to 
the extreme values of the phase space ellipse are to be 
plotted. This will yield a physical representation of what 


the detector will see at the location of interest. 


Zo 





Figure 2-9 2-D Phase Space 


Substituting the system parameters with 1, = .0005, 


Equations M-13 and M-15 become: 


X(t) = .512966, + .739765, +X, (2-3) 


1(t) .073986, + .139406, + .0005 , (2-4) 

where Xo, X(t), lo, 1(t) are in mm; 6, in mrad; (op ue 
+ AP/P. The plot of x(t) vs 1(t) for the points in Fiqume 
2-9 is shown in Figure 2-10. Tabulated values are shown in 
Table 2-2. The center ellipse corresponds to a mono- 
energetic picobunch with divergence in the xz plane. The 
divergence alone causes a Significant increase in the 
picobunch extent. A picobunch of longitudinal extent 1 


micron grows to a longitudinal extent of approximately 14 
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microns at location =e. The right side ellipse represents 
the same picobunch having one-tenth of one-percent greater 


momentum and the left side ellipse represents the picobunch 


With one-tenth of one-percent less momentum. A non- 
monoenergetic picobunch with § = +/- .1% (revised MKIII 


IRFEL parameter) would have a total longitudinal extent of 


= 46 microns at location E. A coherent TR effect would not 
be expected under these conditions. Even in the mono- 
energetic case, the beam divergence would have to be 


greatly reduced to yield a longitudinal extent of <= 3 
microns. The TRANSPORT system equations for detector 


placement midway between locations B and C are: 


Xc(t) = (cosa)x, + (Cesine + Li /coso ee 
+ fo 6 (1-cese)” + by Caricaiieeas (M-5) 
let) = (=Sim ) xo > cet! =e0s oc) -oaaael 
= 06 (Gso ime) 65: (M-7) 


Substituting the known parameters with L, = 163.52 mm these 


equations become: 


x(t) = .98079x, +.204826, + .362784 , (2-5) 


1(t) = -.19509x, - .003756 -.002466 + .0005 . (2-6) 


Zo 


meplet of I(t) vs x(t) for the phase space of Figure 2-9 is 
shown in Figure 2-11. Tabulated values are shown in Table 
>. 

Inspection of Figure 2-11 reveals a shift in the 
orientation of the ellipses with respect to Figure 2-10. An 
11.25 degree rotation of the thin foil about the y-axis will 
Minimize the total longitudinal extent of the original 
picobunch as seen by the thin foil. The effective 
longitudinal extent of the original picobunch is 
approximately 23 microns for constant x(t). Moving the 
detector further upstream to a location 10 mm downstream of 


the first bending magnet the TRANSPORT system equations are: 


pac) = .98079x, + .048309, + .057425, , (2-7) 


it) = —-.19509x, - .003759, - .002465, +.0005 . (2-8) 

Peeolocec Of 1(t) vs x(t) for this detector location is shown 
in Figure 2-12. Tabulated values are shown in Table 2-4. 
Note that Equation 2-6 and Equation 2-8 are identical, thus 
the longitudinal extent of the picobunch 1s constant between 
B and C. However, the effective picobunch extent will 
increase between B and C. Figure 2-12 reveals that a 
momentum spread of one-tenth of one-percent will not cause 


Significant increase in the effective longitudinal extent of 
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thes plcobunci- The non-monoenergetic ellipses appear to 
closely coincide with the central momentum ellipse. 

It is apparent from equations M-7 and M-15 that the 
change in the longitudinal extent of the original picobunch 
is due to the initial horizontal displacement of the 
electron, X59, the divergence of the electron, 6&,, and the 
momentum increment, 5,5; all with respect to the assumed 
central trajectory. A scaled plot is used to illustrate the 
relative effect of each factor contributing to the change in 
the horizontal and longitudinal extent of the picobunch as 
it traverses the system. Substituting the known parameters, 
the TRANSPORT equations to be plotted reduce to the 


equations of Table 2-1. 


TABEE. 21 
EQUATIONS FOR TRAJECTORY PARAMETERS 
AT VARIOUS LOCATIONS 
Detector 10 mm from bending magnet #1: 


x(t) 


-98079xX_5 + .048306, + .057424 
1(t) = -.19509x, - .003756, + lo - .00246% 


Detector halfway between bending magnets #1 and #2: 


x(t) = .98079x, + .204826, + .362786, 


L(t) = -.19509x, -.003756, + lo -.002468, 


Detector at location E: 


x(t) = X95 + .5129665 + .7397656 


1(t) 


.07398 64 + lo + .139405, 
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Figure 2-13 is an unscaled representation of the 
effect of the initial horizontal displacement, x,, on the 
longitudinal extent of the beam. The beam parameters are 
specified by the given 2-D phase _ space. Figure 2-13 
illustrates that the initial horizontal displacement with 
respect to the central trajectory will have an effect on the 
longitudinal extent of the beam between the magnets but will 
have no effect beyond the second magnet at location E. 

Figure 2-14 is a scaled representation of the effect 
of the divergence, 9, on x(t) and 1(t). The effect of the 
divergence term on 1(t) is small over the regions between 
the magnets but increases beyond the second bending magnet. 
Note that the change in the longitudinal extent is measured 
with respect to the beam direction. 

Figure 2-15 illustrates the effect of the momentum 
mmerement, S,, on x(t) and 1(t). 

Comparing Figures 2-14 and 2-15, it appears that the 
divergence and momentum increment have roughly the same 
effect in debunching the beam picobunches as they traverse 
the system. The previous three figures and Figures 2-10, 2- 
11, and 2-12 lead to the physical representation of the 


beam picobunches shown in Figure 2-16. 
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TABLE ®2Z= 
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BEAM PARAEMTERS FOR DETECTOR AT LOCATION E 
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BEAM PARAMETERS FOR DETECTOR HALFWAY 
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BEAM PARAMETERS FOR DETECTOR 10 mm 
DOWNSTREAM OF FIRST BENDING MAGNET 
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III. RESULTS AND CONCLUSIONS 


The axially modulated beam of the MKIII IRFEL was 
represented as disks of charge (picobunches) as shown in 
Figure 1-4. If these picobunches were to retain their axial 
modulation and pass through a thin foil, coherent optical 
transition radiation would be generated. 

The analysis of the TRANSPORT system equations for tk- 
detector locations of Figure 2-5 reveal that significant 
debunching (of more than 3 microns in the axial direction) 
Will occur at these detector locations due to the effects of 
a non-monoenergetic, divergent beam. 

For detector placement at the four-way cross location 
(location E of Figure 2-5) the increase in the picobunch 
extent iS approximately 46 microns (Figure 2-10). For the 
detector location halfway between bending magnets #1 and #2 
(Figure 2-5) the effective picobunch extent is reduced to 23 
mMicrons (Figure 2-11) with the detector rotated 11.25 
degrees about the y-axis. Placing the detector further 
upstream has the overall effect of reducing the effective 
picobunch extent. The axial modulation of the electron beam 
of the MKIII IRFEL can be measured at a detector location 10 
millimeters downstream of the first bending magnet with the 
detector rotated 11.25 degrees about the y-axis. At this 


location there will be a 2 micron increase in the effective 
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peeobunch extent (Figure 2-12). In theory, detector 
placement at this location will enable the desired OTR 
effect to be observed. This location, however, is not 
easily accessible due to structural hardware upstream of the 
bellows of Figure 2-4. In addition, the bellows location 
and waveguide material would have to be modified in order to 
mount a detector of the type illustrated in Figure 2-6. To 
situate a detector at this location it will be necessary to 
make hardware modifications to accomodate the experimental 
arrangement. 

The normalized horizontal and vertical emittance 
parameters for this experimental proposal were 8 mm-mrad and 
4mm-mrad respectively. The momentum spread was one-tenth of 
one percent. A reduction in these parameters will yield a 
reduced picobunch extent at any detector location. 

In addition to the Stanford MKIII IRFEL, there are 
several IRFELS that are potentially suitable for an 
experiment of this nature. Potential candidates for this 
experiment include: 

(1) The Free Electron Laser Driven by the NBS Microtron 


(2) Stanford Superconducting Accelerator/Free Electron 
Laser SCA/FEL 


(3) Lawrence Livermore National Laboratory ATA/TRW 
Wiggler. 


The NBS CW Microtron Driven FEL (currently under 
development) has the design goals of a normalized transverse 


emittance less than 10 mm-mrad and a momentum spread of less 


a2 


than .03% for operation at 0.2-3.0 microns. This momentum 
spread is one order of magnitude better than the MKIII IRFEL 
due to the inherent acceleration process of a microtron. 
The hardware and bending magnet configuration downstream of 
the wiggler exit need to be investigated and modeled in 
order to determine a suitable detector location to measure 
the electron beam modulation. 

The Stanford Superconducting Accelerator/Free Electron 
Laser (SCA/FEL) has a normalized transverse emittance of 5 
mm-mrad (10 mm-mrad with a high brightness injector) and a 
momentum spread of approximately .04% for a fundamental 
wavelength of 3u. [Ref.8] The hardware and bending magnet 
configuration downstream of the wiggler exit also need to be 
investigated and modeled in order to determine suitability 
for a proposed experiment. 

The Advanced Test Accelerator/TRW Wiggler FEL utilizes 
induction linac technology and the laser amplifier FEL 
configuration to amplify a carbon-dioxide laser at 10.6 
microns. This system may be more suitable for detector 
placement due to the absence of optical cavity mirrors and a 
portable beam dump. 

The TRANSPORT equations were generated using a first 
order approximation and an assumed beam spot size of radius 
1 mm at the entrance pole face of the first bending magnet. 
The accuracy of the equations could be improved by utilizing 


the second order fitting capability of TRANSPORT. 
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The system of Figure B-3 was modeled starting at 
location A, assuming a beam spot size of radius 1 mn. A 
more accurate model would incorporate the effects on the 
electron beam envelope within the wiggler. TRANSPORT does 
not have a designated "wiggler" element however, the wiggler 
can be modeled using drift spaces and the "Arbitrary 
Transformation" (type code 14) input option. This model 
should start at the wiggler center where the beam waists Xw 


and yy are known. [Ref. 10) 
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APPENDIX A 


BEAM TRANSPORT THEORY 


1. PHASE SPACE 
The transport of a charged particle beam through a 
system of magnets and drift spaces involves the study of a 
fundamental conservation law known as Liouville's theorem: 
i under the action of forces which can be derived from 
the Hamiltonian, the motion of a group of particles is 
such that the local density of the representative points 
in the phase space remains everywhere constant. [Ref. 11] 
Liouville's theorem advises the reader that he must look at 
the collective group of particles rather than an individual 
particle in order to access the overall nature of beam 
transport. A single particle moving in a three dimensional 
coordinate system is completely specified if its position 
and momentum components are known, namely x, y, Z, Py, Py, 
Po. It is convenient to express this information in a six- 
dimensional space referred to as phase space. A beam 
comprises all individual particles within this phase space, 
neglecting electrostatic repulsion and spin dependent 
effects. vhe latter two effects would require an 
unmanageable phase space of dimension six times the number 
of particles in the beam plus the additional dimensions 
corresponding to the rotational degrees of freedom of the 


individual particles. Any system of conservative forces 


acting on a beam possesses a Hamiltonian. The macroscopic 
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magnetic forces along the beam trajectory of the proposed 
experiment are conservative thus the beam must obey 
Liouville's theorem. The concept of a six-dimensional phase 
space is a valid and useful utility in the study of beam 
transport since space charge effects become weaker as beam 
energy increase. The six-dimensional phase space 
representing all particles in the beam is considered a 
hypervolume that can change shape but must maintain constant 
volume. The phase space analysis of the beam can be further 
simplified by the choice of the coordinate axis. If the 
three components of motion are mutually independent in real 
space, the motion is confined to the three planes (xX, Py), 
(y, Py), and (2, Pg). The six dimensional phase space 
hypervolume has in effect been reduced to three two- 
dimensional phase space areas. A further simplification in 
the analysis is possible by assuming the axial momentum of 
the beam is constant. This is a valid assumption downstream 
of the accelerator sections of the beam trajectory. The 
angular divergence of a particle relative to the beam axis 
is equal to the ratio of the transverse and axial momenta. 
Considering the transverse x-direction the axial divergence 
would be P,/P, or dx/dz. Dimensional phase space (x,P,) is 
now represented as (x,dadx/dz) where dax/dz is the angular 
divergence, a directly observable quantity. Similarly, the 
transverse y-direction phase space representation is 


(y,dy/dz). The longitudinal phase space is represented as 
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(1,dP,/Pz) where 1 denotes the longitudinal beam extent and 


adP,/P, denotes the fractional momentum spread. [Ref.. 11] 


2. EMITTANCE 
The quality of any electron beam is measured in terms 
of its emittance. Consider a Single electron propagating in 


the z-direction as shown in Figure A-1. 





electron having off axis 
component of velocity 









dx/dz=6 


tH electron on beam axis 


Figure A-1 Single Electron Motion in the Horizontal Plane 


Emittance is conveniently expressed in terms of the phase 
space parameters x and dx/dz. Using the small angle 
approximation dx/dz can be represented as 6 Electron beams 
produced by accelerators have elliptical phase- spaces as 


shown in Figure A-2 where: 


Omax = the maximum divergence of an individual electron 


Xmax = the maximum horizontal extent of a beam electron. 
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Figure A-2 Elliptical Phase Space 
In accordance with Liouville's theorem, phase space is 
conserved. The electrons are within the perimeter of the 
ellipse. Emittance is mathematically defined: 

€E = (area of phase space)/r . (A-1) 
Normalized emittance is defined: 

En = Oe, (A-2) 

where y is the ratio of the particle energy to the rest mass 


energy and B= 1. The area of the elliptical phase space is 


1tXmaxé8max therefore: 
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En = Xmaxomax. 


For example: The Stanford MKIII FEL has the following 


normalized parameters [Ref. 12]: 


horizontal phase space area = 8 tmm-mrad , 


vertical phase space area = 47 mm-mrad . 


therefore assuming E = 40MeV: 


-10 mm-mrad , 


Enh = 8 mm-mrad; ep, 


-0O5 mm-mrad . 


Envy = 4 mm-mrad; ey 


Assuming elliptical phase spaces, we can assign approximate 
values for Xmax and Ymax and the corresponding values for 
6émax and ¢max. 

Letting Xmax = Ymax = 1mm and energy equal to 40MeV with 
oe 1 the corresponding values of max and max are 


determined to be: 
XmaxX = 1 mm; @6max = .10 mrad , 


Ymax = 1 mm; gmax = .05 mrad . 
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3. BEAM WAISTS 

In accordance with electron beam transport design [Ref. 
10], the electron beam is normally focussed at the center of 
the wiggler of the FEL. This yields a beam waist in the 


center of the wiggler as shown in Figure A-3. 


Figure A-3 Beam Focus at Center of Wiggler 


where 


Xy = horizontal beam waist 
Yw = vertical beam waist 
Ly = length of wiggler (1 meter) 


The horizontal beam waist size is related to the length of 


the wiggler and the horizontal emittance by the relation: 


Xyw = Lwe h/ 2 (A-3) 
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The vertical beam waist size is related to the vertical 
emittance and the intrinsic betafunction, Bw, by the 


relation: 





Yw = Vey 8y (A-4) 


where 


Bw = electron radius of curvature at 
center of wiggler. 
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APPENDIX B 
TRANSPORT-A COMPUTER PROGRAM FOR DESIGNING 
CHARGED PARTICLE BEAM SYSTEMS 
1. INTRODUCTION TO TRANSPORT 

The study of static-magnetic beam transport systems led 
to the development of the original first-order TRANSPORT 
computer program. TRANSPORT was written in BALGOL at 
Stanford Linear Accelerator Center (SLAC) in 1963. It has 
been refined in the succeeding years by contributions from 
various laboratories throughout the world. Subsequent 
improvements include the translation into FORTRAN and a 
second-order fitting capability. The version of TRANSPORT 
used in this thesis is the PC version dated May, 1985 from 
D. Carey, Fermi National Accelerator Laboratory (FNAL), 
compiled with the IBM version of professional fortran by 
Bvan—-McFarland. [Ref. 13] 

TRANSPORT utilizes matrix theory methods to step 
sequentially through a beam line system consisting of 
magnets, drift spaces, and special magnet configurations. 
The beam line system and the input beam parameters are 
specified by the user and an output beam matrix and overall 
system transfer matrix are generated. Each individual step 
through the system 1s concidered an element. Each element 
is given by a sequence of items separated by spaces and 


terminated by a semicolon. These items are a type code 
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number, a vary field, the physical parameter and optional 
labels. The type code number identifies the element (i.e., 
magnet, drift space, pole face rotation). A vary field 
allows the user to adjust the physical parameters of an 
element if the problem involves fitting. Physical 
parameters are quantities which describe the physical 
element. Optional labels assist in identifying beam 
elements and parameters in the output. A summary of 
TRANSPORT type codes is shown on page 77. {[Ref. 13] 

The initial objective in writing the input program will 
be to specify the beam line, establishing the central 
Erajyeceonuys The system to model by TRANSPORT is shown in 
Figure B-l. This system consists of two bending magnets and 
five drift lengths where 

Lo = initial drift length , 

L = effective pathlength through bending magnet , 
yi = eee emg Ly 

Lo = drift length 2 , 


a = bending angle . 


ae SAMPLE PROGRAM DEVELOPMENT [Ref. 13] 

Figure B-2 is a sample input program that establishes 
the beam line of Figure B-1, assuming the initial drift 
space, Lo equals Zero. Refer to Figure B-2 for the 


renarnadeseot Ehis section. 
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¢— bending magnets 





Figure B-1l System Configuration 


ot 


Cr ‘STANFORD MKITI IRFEL: REVISED [O07 te73i 


0 
eZ) (CARTESIAN COORDINATE SYSTEM; x,y,z) 
1:3), “COOR- Peo, : 
(2) (GET UNITS/scaled units) 
Loi *xymm' joy * mm’ 0 eet ; 
bey: "EMeV’ ie "MeV' 07001 
Se eerie G. ’ mim’ C200T 
tou ’leng’ oe a Seni 0.41: Je 
(4) (BEAM PARAMETERS AT FIRST BENDING MAGNET POLE-FACE) 


(HORMALIZED BEAM EMITTANCE: Enh=8mmn-mrad, Env=4mm-mrad) 
20 ’ BEAN’ 1.0 Orr 8 0.05 0005 0-0 40 01n 
(monoenergetic beam) 


(5) Lo. 52 520 ; (half-gap). 
(defanlt value of K1=.5) 
(6) (FIRST BENDING MAGNET ) 
ree Gir : 
Ne. "BND ’ Sonat 6.84175 O: ‘ 
ae PL eZe ; 
ie re * (PREWT PRAM MATRIX) 
GS 4. > (PRINT XFORM WATRIX) 
Gr) (DRIFT GPACE #1) 
3. "DETL S200 
ese 122 ’ 
tele 4. : 
(8) (SECOND RENDING MAGHEET) 
20. 180. *LEET” (std coord rot w/11.25 bend) 
a Wee Ol ten 
4, "BND? ’ 33-347 6.63175 0. : 
ae QO. : 
at), =), : 
on 1 ; 
Pas te ; 
(9) (DRIFT SPACE #2) 
3S “eae 5 80 Ss; 
ee po : 
1. 4. ; 


(70) SENTINEL  ; 


Figure B-2 Sample Input Program 
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The individual elements and there associated type codes will 
be explained in the succeeding paragraphs. 

(1) ‘STANFORD MKIII FEL: revised 10/16/87' is the 
arbitrary title of the sample input program. The zero 
immediately below the title indicates the start of a new 
problem. 

(2) A cartesian coordinate system was chosen for the 
proposed experiment design problem. The coordinate system 
is established using the type code 13. The statement [13. 
12. ;] will yield a beam line reference trajectory that will 
start at the origin and proceed along the positive z-axis. 
The y-axis will point up and the x-axis will point to the 
wert. Statements enclosed by ( ) are comment statements to 
assist in following the code. 

(3) The scaling of units is accomplished via type code 
nes All lengths in Figure B-1 are to be specified in 
millimeters and energy is to be specified in MeV. TRANSPORT 
has the standard set of units as shown on page 78. 
Conversion factors are applied to code digits 1, 11, 8, and 
5 to establish an appropriate unit system. A summary of 
TRANSPORT type codes is found on page 77. 

(4) The input beam, denoted by the type code 1, is 
specified in terms of its phase space and average momentum 
parameters. The phase space input is given in terms of the 
Semi-axis of a six-dimensional beam ellipsoid representing 


the phase space variables x, 90, y, », 1, and & where 
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x = one-half the horizontal beam extent (1.0mm) , 


0 = one-half the horizontal beam divergence (0.lmrad) , 
y = one-half the vertical beam extent (1.0mm) , 

d= one-half the verticle beam divergence (.05mrad) , 

1 = one half the longitudinal beam extent (.0005mm) , 

6 = one-half the momentum spread (0.0%) , 

P= average momentum of the beam (40MeV/c). 


NOTE: The above units were scaled using type code 15. 

(5) and (6) Bending magnet parameters are specified in 
terms of the wedge magnet, type code 4; pole-face rotation, 
type code 2; and a gap setting, type code 16. A wedge 
bending magnet implies that the central trajectory of the 
beam enters and exitS perpendicularly to the pole face 
boundaries. A wedge magnet can be modified into a square 
edge magnet using the pole face rotation element on either 
the entrance or exit face of the magnet. The field 
boundaries for bending magnets are shown on page 73. ie 
modify a wedge bending magnet into a square for a bend to 
the “right, B, is set to zero and 85 is set to 11.25 
degrees. The pole face rotation type code 2. must 
immediately precede or immediately follow the bending magnet 
(type code 4.0). There are four first order parameters to 
be specified for the wedge magnet; the type code, the 
effective length of the central trajectory, the central 
field strength and the field gradient. The bending magnet 


equations and parameters for the system are listed below: 
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L = 9oa where 


L = effective pathlength through magnet , 

09 = bending radius; (195.30 mm) , 

a = bend angle; (11.25 degrees) , 

L = (195.30mm) (11.25deg) ( rad/1i80deg) = 38.347 mm 


(NOTE: To Simulate a gradual variation of the fields at 
the ends of the bending magnets a fringing field 
element is provided.) 


BO) §=7337.3596(P/0,) where 
B(O) = central field strength , 
P = central momentum; (40 MeV/c) , 
B(0O) = 33.356(40MeV/c) /195.30mm = 6.83175 kG , 
n = magnetic field gradient (n = O assumed). 


Code [{13. 1. ;] yields the printing of the current beam 
feo matrix. Code [13. 4. ;] yields the printing of the 
current transformation matrix, Rl. These codes will be 
addressed in the matrix formulation of TRANSPORT. 

(7) and (9) A drift length (type code 3.0) is a field 
free region through which the beam passes. It is denoted by 


the type code and the effective drift length: 


woe ert PTength 2 (327.03 mm) , 
Lo = drift length 2 (95.30 mm) 
(8) The second bending magnet differs from the first 


in that 8, is set to 11.25 degrees and & is set to zero 
degrees. In addition, a coordinate rotation element (type 
code 20.) is required to rotate the transverse coordinates x 


and y through an angle about the z-axis. A bend to the left 


55 


(looking in the direction of beam travel) 1s accomplished by 
rotating the x,y coordinates by 180 degrees. 
(10) The word "sentinel" denotes the end of the input 


program. 


3. MATRIX FORMULATION OF TRANSPORT (Ref. 13] 
A. Phase Space Vector 
At any specified position in the system an arbitrary 
charged particle is represented by a vector xX, whose 
components are the positions, angles, and momentum of the 
particle with repect to the central trajectory. The X 


vector is represented: 


x 
9 
yi 
X = o 
uy 
| 6 
where 
x = the horizontal displacement of the arbitrary ray 


with respect to the assumed central trajectory , 


6 = the angle this ray makes in the horizontal plane 
with respect to the assumed central trajectory , 


y = the vertical displacement of the ray with respect to 
the assumed central trajectory , 


¢ = the vertical angle of the ray with respect to the 
central trajectory , 


1 = the path length difference between the arbitrary ray 
and the central trajectory , 


= AP/P; the fractional momentum deviation of the ray 
from the assumed central trajectory. 


O> 
| 
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B. The Transfer Matrix R 
The passage of a charged particle through an element 
(magnet, drift space, rotation, etc.) may be represented by 


Equation B-1: 


X(1) = RX(0) , (B-1) 


where 


X(0) is the initial coordinate vector, 
X(1) is the final coordinate vector, 


R is the transformation matrix of the element. 


In a system of several magnets, pole face rotations, and 


drift spaces, R is replaced by the product matrix R(t): 


Ree Rim) ooeR(S)R(2)R(L) « (ae) 


R(t) is the product of the individual element matrices of 


the system. Equation B-3 represents the matrix equation for 


a system: 
6 
Xi(t) = ) Rij%j 09) , (B-3) 
i 
where 
a 
Xo = 6G, 


a 


~ 
ee) 
ll 
M< 


This product matrix is automatically calculated by the 
program and is called TRANSFORM1. TRANSFORM1 can be printed 
when desired using the statement [13. 4. ;] (e.g., paragraph 
(6). =: For static magnetic systems possessing midplane 
symmetry, the six simultaneous linear equations renresented 


by Equation B-3 may be expanded in matrix form as [Ref. 14}: 


x(t) epoca O O Ri¢ “_ 
e(t) R91 Roo 0 0 0 R 6 a 
y(t)} | 0 0 R33 R3q 0 0 ve 
o(t) O 0 Rigs Rage 0 b6 
Ge) R51 R52 0 0 1 Rs¢ 1 
§(t) O O O 0 O a 50 | 
L J 
The transformation is from an initial position T = 0 
to a final position T = t measured along the assumed central 
trajectory. The final coordinate vector X(t) is the matrix 


product of the individual coordinate vector and the overall 
system transformation matrix. The individual Ris elements 
of the matrix are dimensionally consistent with the input 
units specified in the transport code. The individual 


transformation matrices needed to represent the system of 
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Figure B-1 are shown on pages 74-76. It is desired to 
determine and plot x(t) vs 1(t) for the overall system. 
Since the system 1S symmetric, x(t) and 1(t) have no 
dependence upon yo, or %, thus the general form of the 


tranformation matrix is reduced to: 


x(t) Rii Rig 0 Ras | Xo 
0(t) Ro1 Ro2 0 R26 86 
1(t) R51 R52 1 R56 lo 


O 
be 
— 
Or 
0 
t 


s(t) | | 0 


The wedge bending magnet is reduced to: 


cos a OoSin a 0 Pg (1-cose ) | 
| 

-hsing cos a 0 sin a 

-Sin g Bion L=Cos G) il = Oay( ese ue 





— 


The pole face rotation matrix is reduced to: 


1 0 O O 
(tam )/po6 sli 0 0 
0 0 dl 0 
bh i) O O i 
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The drift space matrix is reduced to: 


a L 0 0 
0 Al 0 0 
0 0 x 0 
6) 0 0 1 |. 


It is desired to determine the X matrix at various locations 
along the central trajectory. Figure B-3 denotes the system 


locations to be evaluated: 


¢— bending magnets 





Figure B-3 System Configuration Evaluation Points 
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The matrix representation of Xp(t) is the product of the 
pole face rotation matrix, wedge bending magnet matrix, and 


the initial coordinate vector: 





Bee ( C) Ih Cl 0) 10, Cos a p551n0 0 29(1-cosa) || xo 
6B(t) calle, p54 70 0 ~hsina cosa 0 sind | a 
a(t ) 0 Oe. Fy 70 SECM L-COs ca) a -poCe-8ina)) 19 
S(t) 0 0 0 1 0 or 1 bo. | 
(pole face rotation) (wedge bending magnet) 


The matrix product reduces to: 








Xp (t) cos Oo O9S1in a 0 PelGlS cosa) [ a | 
6 p(t) 0 iP eos 0 tan o a5 
ie ( c) 7 -sing = 09 (1-COSs q) 1 aco or Sama) i | 
6 R(t) 0 0 0 il oz | 


Biewmatrix representation of X(t) is Xp(t) multiplied by 


the first drift space matrix: 


oa Fi 


xc (t) ee Ono cosa sina 0 [% (1-cos a) xe 

9 a(t) O al o_O 0 1/fcosa 0 tan a oe 

le(t) y | O O 1 0 aoe cost) 1 “= 5. (o-Sin o lo 
| sect) | | 0 0 0 1 | 0 0 0 i | oF il 
(first drift) ees cs : 


Gar 


The matrix product reduces to: 


| xc (tC) cosa P)>sSinatL,/cosa 0 Po(1l-cosa)+L,tan a | Xo 

| 

| 8e(t) 0 1/cosa 0 tan a Oe 

: ie Ce) -sina => (1-eesa) 1 -09(G=-sina@) lo | 
| S(t) | 0 0 0 ao ; So | 


The matrix representation of Xp(t) includes a coordinate 
rotation which will in effect bend the beam to the left. [In 
matrix form this is represented by changing the bending 
radius and bending angle to its negative value (1.e., 0g -~ 
“9, & —a)- The wedge bending matrix for a bend to the left 


is modified as shown 


cos a PoSin a 0 05 (Secs) 
-hsina cosa 0 -sina 
sing 06 (1-cos q) di —po (cqs lhe 


0 0 O als 


The pole face rotation matrix retains its identical form 
because tan(-q) = -tan(q). Multiplying Xc(t) by the pole 
face rotation matrix and the wedge bending magnet matrix 


yields: 
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Xp (Cc) i" 2 og tan ot Ly /cos? a 0 z1 | | Xo | 


| 
8p (t) 7 0 1 0 oa | 86 
TD ©) 0 Te } Ze lo | 
Sptt) | | o : 0 0 il 86 | 
u ‘ _ 
where: 


Zl ==-205 + 2P,/cosa + L,sino/cos4 a ; 
Z2 = 21 , 
(oe =2 00 + 205cand + Li tan? a . 


The matrix representation of Xp(t) is Xp(t) multiplied by 


the second drift space matrix: 





[xp (t) 1 Z4 0 Z5 Xo 
a(t) 0 a 0 0 c= 
le (t) 0 Z6 1 Ze | a 

| ect) | 0 0 0 dl . De 
where: 


24 = 20,tana + Ly/cos? a + Lo , 
eee eet 2 05/ COS Or Lisino/cos¢a , 
Ga 25 3 

ae ented eee aly Ot bytance o . 
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In equation form, the matrix products are represented as: 


Maat) (ceosa)x,4 + (positmame. a 0 (1-cos a) So (M-1) 
Gp (t) (1/cosa) 66 + (tana) & (Me 2 
tate) (-sing)X5 - @(1-cOSa) 866 + lo - O(a — Sind) SO 
6p (t) 50 (M-4) 
Xc(t) (cosa)X5 + (p5Si2na + Li/cosarg, eit —ecos (Ms 
+ Lytana)] & 
Ga (t) (1/CoS 0) 6 + (tan®) 5, (M-6) 
la(t) (-sind) xg - %(1-cosa) & + lo - P6(%sin® 6 (M-7) 
CeCe) O65 (M-8) 
Sock) X_(t) + Ly 6g(t) (M-5A) 
fo (t) Og (t) (M-6A) 
1e(t) ee) (M-7A) 
6c(t) =6 p(t) (M-8A) 


Xp(t) = xo + (2 pecan as L,/cos¢  & + (-2% + 20%/cos (M-9) 
+ Lysin o/ cos a) & 


Die = clo (M-10) 


dy CE) 2S 2 Ge 2G COs a bean of cos? Q8o + Be (M—11) 
+ (-2p50 + 2potana + Lytan7a)6, 
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Baye) = (2/cosa)xc(t) + oot c(t) — 5 (L-cos o) 6c(T) 


(M-9A) 
P(t) = (cosa) 6e(t) - (sina) 6¢(t) (M-10A) 
Ip(t) = (tana)xce(t) + (eee Osajogm (eG) + le (M-11A) 

Src, — Salat sens), 
Sp(t) = %e(t) (M-12A) 
es tec cot)s, >” 
E(t) = %& (M-14) 
lp(t) = (-2p9 + 2@/cosa+ Lysim/cos? 6, + 1, (M-15) 
P(=2050 eZ eocane +L, tan2 a) 8, 

Sa(t) = % (M-16) 
Pee) = Xp(t) + Lo h(t) (M-13A) 
On (t) = (tt) (M-14A) 
Ip(t) = 1p(t) (M-15A) 
8B(t) = Sp(t) (M-16A) 


C. Beam Matrix 

The phase ellipse beam (sigma) matrix is used to 
inform the user of the beam parameters at any location in 
the system. The beam matrix will represent the bundle of 
particles or rays constituting the phase space of the input 
‘beam at desired locations in the system. Particles of the 
beam are assumed to lie within the boundaries of the 
ellipsoid with each point within the ellipsoid representing 
a possible ray. The equation of an n-dimensional ellipsoid 


may be written in matrix form as: 
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X(0) 4g (0) == xe) eee (B-4) 


where: 


x(0)? = transpose of coordinate vector X(0), 
o6(0) = a real positive definite, symmetric matrix. 
As a particle passes through the system, it 


undergoes the matrix transformation of Equation B-1. 
Combining this transformation with the equation of the 
initial ellipsoid and using RR71 = I 

x(0)?criRf-1) o(0)~1(R71tR)x(0) = 1. (B-5) 
Equation B-5 can be rewritten as: 


[RX(0)}]?{Ro(0)R?)~1(RK(0)] =1 . (B-6) 


The equation of the new ellipsoid after the transformation 


becomes: 


X(1)To(1)71x(1) = 1, (B-7) 


where: 


o(1) = Ro(o)R? . (B-8) 
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Thus the beam may be represented at any point in the system 
by Equation B-8. Note: o(0) is the input beam phase space 
specified by the user. The beam matrix is printed using the 
statement [{13. 1. 7]. The projection of the semi-axis of 
the ellipsoid upon each of its Six coordinate axes is 
printed in a vertical array, and the correlations among 
these components are printed in a triangular array. The 
phase ellipse beam matrix in the printed output form is 


shown in Figure B-4. 





x Y3(11) CM 

8 ¥O(22) MR reel) 

y 70 (33) CM r(31) r(32) 

> 7o(44) MR r(41) r(42) r(43) 

£ 7o(55) CM r(51) r(52) EC53) r(54) 

5 4c (66) PC r(61) r(62) r(63) r(64) r(65) 
where: TCs) ep OE : 


[o(iiye(ji)] 3 


As a result of the fact that the O matrix is positive definite, the 


r(ij) satisfy the relation 


Cees ls 


Figure B-4 Beam matrix 
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(Note: units will be labeled as specified by type code 15.) 
The individual sigma elements are defined: 
Yo (11) = Xpax = the maximum half-width of the beam 


envelope in the x(bend)-plane at the 
point of Ehe pEIMntoue- 


Yo (22) = Omax = the maximum half-angular divergence of 
the beam envelope in the x(bend) plane, 

Yo (33) = Ymax = the maximum half-height of the beam 
envelope, 

vo (44) = dmax = the maximum half-angular divergence of 
the beam envelope in the y(non-bend) 
plane, 

Vo (55) = lmax = one-half the longitudinal extent of the 
bunch of particles, 

Vo (66) = 64 = the half-width 1/2(4P/P) of the momentum 


interval being transmitted by the system. 


Note: the units appearing next to the Yo (ii) in jeme 
TRANSPORT printout are the units chosen for coordinates x, 
GQ. Vio; 19 and. 6 respectively. 


An example of the (x,98) plane ellipse is illustrated 
in Figure B-5. 


Consider a two dimensional (x, 6) plane projection 


of the general 6-dimensional ellipsoid. Let: 


L Cia O21 
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Figure B-5 A Two Dimensional Beam Phase Ellipse 


be a real, positive definite, symmetric matrix. The inverse 


of the O0-matrix exists and is defined by: 
022 mo 
mec ean 
where €* is the determinant of . 
The two-dimensional coordinate vector and its 


transpose are represented by: 


X = eo 76) = 
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The expansion of the matrix equation 
xP oclxy = 1 
is the equation of an ellipse, therefore 
DO 55x* - 206 51x9 + 01, # = €4 = deto. 


The area of the ellipse is given by: 


2 
A =7 (deto ) = TXmax Sint = ™Xint’max = Te- 


The correlation between x and 8 (the orientation of the 
ellipse) depends on the off diagonal term O5,- This 
correlation is defined as: 


=o 


eee 


SDA Seekeeoe 
The correlation, r, measures the tilt of the ellipse and the 
intersection of the ellipse with the coordinate axis. The 
beam ellipse in the (x, °°) plane will be valuable in the 
analysis of the beam divergence at various points in the 
system and also in the fact that it reveals the change in 


the longitudinal extent of the beam via element Vo (55). 
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4. 


INSTALLATION OF TRANSPORT ON THE IBM PC [Ref. 15] 


The following requirements and procedures are necessary 


for the loading and running of TRANSPORT: 


Requirements: 


config.sys file must have files >= 15; buffers >= 15. 
disk set A containing .EXE and .BAT files in DOS/BACKUP 
640K memory 


hard disk; 10M-byte or larger 


Procedures: 


move to the C directory: C: [ENTER] 
create a TEMPO subdirectory: MD\TEMPO[ENTER] 
move to the TEMPO subdirectory:  CD\TEMPO[ENTER}] 


RESTORE .EXE and .BAT files: RESTORE A: C: [ENTER] 
(insert disk #1 and #2 as directed) 


move to the C directory: CD\C: [ENTER] 

create a TRANS subdirectory: MD\TRANS[ENTER] 

move to the TEMPO subdirectory: CD\TEMPO[ENTER] 

copy files from TEMPO to TRANS: copy *.* \TRANS[ENTER] 
delete files from TEMPO subdirectory: del *.*[ENTER] 
move to the C directory: CD\C: [ENTER] 


remove TEMPO subdirectory: RD\TEMPO[ENTER] 


Limitations: 


the PC version of TRANSPORT will allow for 1000 elements 
and 4500 data values 


accuracy is four significant figures 


7 a 


To° run TRANSPORT: 


- TRANSGO.BAT 1s the command file to run TRANS.EXE 


- type the following: 


TRANSGO FILES FILEZ Pees 


FILE1L: input data file name (users file name) 

FILE2: output print file name (users file name) 

FILE3: output summary file name (users file name) 
- example: 


TRANSGO MK3FELA MK3FELA.OUT MK3FELA.JNK[ ENTER] 


* Note: TRANSPORT is currently loaded in the IBM PCs in the 
Physics Library and Physics Department Office under 
the subdirectory C:\TRANS. X. Maruyama currently 
holds the TRANSPORT software on 5-1/4 disks. 
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Table 1: Summary of TRANSPORT type codes 


PHYS TYTE 3rd ath Sth 6th 7th ath Sth 


& (percent ) 


r.m.s. ADDITION TO 


BEAM ENVELOPE Aé{percent) 


ANGI F .OF 


VOLE FACZ ROTATION RITALION (degrees) 


IRI FT LENCTH (metres) 


FIELD GRADIENT 


BENDING MAGNET LENGTH (metres) FIELD (kG) (aeeetue) 


QUADRUPOLE LENCTH (metres) FIELD (kG) HALF-APERTURE (on) 


TRANSPORM 1 UPDATE 0.0 1.0 

TRANSIUM 2 UPIATR 0.0 2.0 

BLAM CUNTROID SUIT SIFT (x) (on) MET (8) (mr) SMET Cy) (om) SHIFT (0}(mr) SuLET (£) (cn) SHINT (percent) 
ALIGNMENT TOLERANCE DISPLAG2ENT (x)(cm}| ROTATION (8) (mr) DISPLACDENT (y)(oni] ROTATION (¢)(ar) | DISPLACEMENT (2){cm) ROTATION (a) (mr) 


NUMBER OF 
REPEAT CONTROL REPEATS 
UESIRUD VALUE OF 
(1,3) MATRIX 
HMENIS 


ACUURACY OF 
iT 


FITTING QUNSTRAINTS 1 





Note: 1 ls used for fltting a beam (oc) matrix element, “1 is used for fitting an Rl matrix element. 
- {1 * 20) is used for fitting an R2 matrix element. 


E {energy gain) @ (phase lag) 
ACCELERATOR 11.0 LENGTH (metres) (GeV) (degrees) (WAVELENCTHN (cn) 
Fieatéa eltipsey | 22-9 THE FIFTEEN CORRELATIONS AMONG THE SIX ELEMENTS (shis entry mist be preceded by a type code 1.6 entry.) 
INFUT/OUTPUT 13.0 CONTROL 
OPTIONS 5 COLE NUMBER 
ARBITRARY R MATRIX | 14. vvvwwv0j R(J,1) R(J,2) R(J,3) R({J,4) R(J,$) R{J,6) 2 
UNTTS CONTROL : , , 
eer 15.0 cnne SYMBOL eee 
Dimensions) at (af required) 
la 16.0v 1.0 e(1l) « a(t)’ o,in units of transverse length (an) 
Parcan PARTICLES 16.0 3.0 M/m (dimensionless)/ m © mass of alectron 
HALF-APERTURE OF 
BENDING MUNET 16.0 4.0 w/2 {ca) 


IN x-PLANE 

HALF-APIRTURE OF 

EINDING MUXET 16.0 $.0 g/2 (on) 
IN y-PLANE (gap) 

LENGTH OF SYSTEM 16.0 6.0 L (metres) 


FRINCE FIFLD GOR- 





| 
! 

RECTION COEFFICIENT 16.0 7.0 Ky (dimenslonless) | | 

FRINGE FIFILD COR- | 

| CURVATURE OF . | 

ENTRANCE FACE OF 16.0v 12.0 {I/Ry) (L/metres) 

RENDING MAGNET ' 

CUPVATURE OF 

EXIT FACE OF 16.0v 13.0 {I/Ry) (l/metres) 

BENDING MAGNET 

FOCAL PLANE, 16.0 15.0 Angle of focal plane rotation (degrees). 

RUILATION : See type cate 16.0 for details, 

(INITIAL ALAM LINE 

2 -COURDINATE 16.0¥V 16.0 He 

INITIAL BEAM LINE 

INITIAL BEAM LINE 

INITIAL ALAM LINE 

MORIZONTAL ANGLE | '8-O¥ | 19.0 0 

INITIAL BEAM LINE 

SECONT-ORDER 17.0 

CALCULAT TORS ; 

SE XTUPOLE 18. 0v LENGTH (metres) FIELD (kG) HALF-APERTURE {cm) : ' 

SOLENOID 19,.vv LFINCTHE {metres} FICLD (kG) | 

ANGLE OF 
BEAM ROTATION 20.¥ ROTATIUN (degrees) 
STRAY FIELD 21.0 See lster section of report, 








os " a? os a 


Note: The v's following the type codes indicate the parameters which may be varied, See section under type code 10.0 for a detailed explanation of 
Vary Codes. The units are standard TRANSPORT units (as shown) unless changed via type code 15.0 entries. 


qT 


The various units that may be changed are: 


Code Quantity Standard Symbols used in 
Digit TRANSPORT Unit SLAC-75 
1.0 horizontal and vertical trans- cm xy 


verse dimensions, magnet 
apertures, and misalignment 


displacements 


2.0 horizontal and vertical angles mr 6,¢ 
and misalignment rotation 


angles 


* 
3.0 vertical beam extent (only) ) cm y 


and bending magnet gap height 


* 
4,0 vertical beam divergence mr d 
(only ) 
5.0 pulsed beam length and wave cm f 
length in accelerator 
6.0 momentum spread percent (PC) 6 
7.0 bend, pole face rotation, and degrees (DEG) 


coordinate layout angles 


8.0 length (longitudinal) of ele- metres (M) t 
ments, layout coordinates, 


and bending magnet pole face 


curvatures 
9.0 magnetic fields kG B 
10.0 mass electron mass m 
11,0 momentum and GeV/c p(0) 
energy gain in accelerator GeV AE 
section 


*) These codes should not be used if the coordinate rotation (20. 0) type code is 
used anywhere in the system. 
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